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We derive the Bogoliubov+U formalism to study the thermodynamical properties of the Bose- 
Hubbard model. The framework can be viewed as the zero-frequency limit of bosonic dynamical 
mean-held theory (B-DMFT), but equally well as an extension of the mean-held decoupling approx¬ 
imation in which pair creation and annihilation of depleted particles is taken into account. The 
self-energy on the impurity site is treated variationally, minimizing the grand potential. The theory 
containing just three parameters that are determined self-consistently reproduces the T — 0 phase 
diagrams of the three-dimensional and two-dimensional Bose-Hubbard model with an accuracy of 
1% or better. The superhuid to normal transition at hnite temperature is also reproduced well and 
only slightly less accurately than in B-DMFT. 


PACS numbers: 71.10.Fd, 02.70.Ss, 05.30.Jp 

I. INTRODUCTION 

The properties of cold atomic gases trapped in an opti¬ 
cal lattice can be tuned and controlled very precisely, pro¬ 
viding a powerful tool for the simulation of the iconical 
low-energy effective Hamiltonians of condensed-matter 
models [I]. Dramatic experimental progress in this held, 
such as the observation of the Mott insulator to super¬ 
huid phase transition in the Bose-Hubbard model [1. or 
the recent realization of the Hofstadter model |2j, have 
galvanized the condensed matter community. 

The accuracy of cold atom experiments has laid bare 
the need for advanced numerical methods in order to 
tackle these correlated many-body systems quantita¬ 
tively. In one dimension the density matrix renormaliza¬ 
tion group (DMRG) [3H6| works very well, while in higher 
dimensions the numerical complexity represents a prob¬ 
lem. Very successful also in higher dimensions have been 
path integral quantum Monte Carlo (QMC) simulations 
with worm-type updates [7] leading to identical results 
for the Bose-Hubbard model as observed in experiment 
|8I9| . Despite all its impressive successes for bosonic sys¬ 
tems the worm-algorithm suffers from a prohibitive sign 
problem when cold atoms are subject to (artiflcal) gauge 
fields. In such cases no algorithm is known that works 
and one is forced to resort to approximations. This has 
been a main driving force for the development of bosonic 
dynamical mean-field theory [lQIllj (B-DMFT). 

In B-DMFT, as in standard mean-field theory, the 
many-body system is projected onto a single impurity, 
keeping only the local propagators. This provides us with 
a model in which the self-energy and the local propaga¬ 
tors can be determined self-consistently by solving an 
effective impurity action and a self-consistency condition 
iteratively. At the moment it has only been formulated 
for single-site impurities, but the ultimate goal is to for¬ 
mulate it for small clusters that can also deal with larger 
unit cells. It is known that B-DMFT provides excellent 
agreement [of the order of 1% in three dimensions (3D)] 
with experimental and QMC data m for the standard 


Bose-Hubbard model and improves remarkably on static 
mean-field theory. B-DMFT is hence a promising can¬ 
didate to deal with more complicated dispersions: The 
hope is that it deals with the local interactions as accu¬ 
rately as in the standard Bose-Hubbard model whereas a 
small cluster could treat the full dispersion. Furthermore, 
one would expect that a cluster method goes beyond real- 
space DMFT for multi-component systems and systems 
with long-range interactions and/or dispersions (12113) . 

However, B-DMFT is numerically rather complex due 
to the imaginary-time dependency of the hybridization 
term. At finite temperature the impurity problem has 
to be solved by continuous time Monte Carlo methods 
mm, where, due to the difference in sign between the 
normal and the anomalous Green’s function, a sign prob¬ 
lem arises in the symmetry broken phase. 

In this paper, we filter out the ingredients of B-DMFT 
that are indispensable for its accuracy and arrive at a 
simpler formalism. This is the Bogoliubov+U theory 
(B+U), which makes use of a simplified effective impurity 
Hamiltonian, similar to the action of extended mean field 
theory, which was recently developed in the high-energy 
community mm but differs conceptually from our for¬ 
malism. B+U has a negligible computational cost and is 
not prone to numerical instabilities. The premise of our 
theory is that the Bose-Hubbard model can be fully char¬ 
acterized at zero temperature by the three parameters </>, 
En, and E 12 (the condensate order parameter, and the 
normal and the anomalous self-energy, respectively) if 
the self-energy is treated as a variational parameter, pro¬ 
viding a far better approximation to finite-temperature 
properties than simple mean-field theory. B+U can be 
seen as a simplified B-DMFT where only a single Mat- 
subara frequency is kept (and is hence conserving). It 
is different from the variational cluster approximation 
(VCA) by also considering non zero values of pair cre¬ 
ation and annihilation of depleted particles PS- It is also 
the simplest accurate extension of the weakly interacting 
Bose gas theory D2I to lattice systems with a superfluid 
to Mott insulator transition. It further provides a very 
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natural framework compared to the collective quantum 
field theory developed by Kleinert et al. [18] and the 
two collective fields proposed by Cooper et al. USM], 
and behaves quantitatively much better. It has a similar 
functional degree of freedom as the projector technique 
introduced by Trefzger and Sengupta 22] for finite lat¬ 
tices. Among other methods which are used to treat the 
Bose-Hubbard model are also the process chain approach 
[23124] and the quantum phase field U( 1) rotor field [[25] . 

The paper is organized as follows. In Sec. [n]the B+U 
formalism is derived for the Bose-Hubbard model in equi¬ 
librium, while in Sec. |III| the details of the variational 
calculation of the optimal self-energy are shown. We 
furthermore summarize the full self-consistent scheme of 
B+U and show how thermodynamic quantities can be 
calculated from it in Sec. |IV[ while some simple limits of 
B+U are explained in Sec. |V| In Sec. [Vl]we present the 
results at zero and finite temperature comparing them 
with QMC and B-DMFT. Finally, in Sec. |VII| we con¬ 
clude and present a short outlook about future applica¬ 
tions of the B+U formalism. 


II. SOLVER AND SELF-CONSISTENCY 
CONDITION 


In this section we derive the B+U formalism for the 
Bose-Hubbard model in equilibrium. In order to derive 
an effective Hamiltonian, we start from the full Bose- 
Hubbard Hamiltonian, 

H bh = - J £ b\bj + (!) 

(ij) i i 


where b\ is a bosonic single-particle creation operator on 
lattice-site i, n i is the particle-number operator, J de¬ 
notes the tunneling amplitude, U the on-site interaction, 
fi the chemical potential, and (i,j) means that we sum 
over nearest neighbors. In order to determine the ther¬ 
modynamic properties of the system, we have to compute 
the condensate density and the connected Green’s func¬ 
tion, defined respectively as 

<fi=(b), (2) 

G^'(t) = - ($,(r)$]( 0 )) + U\ ( 3 ) 


with Nambu notation 4q(r) = 


bi{r) 

b l( T ) 


The possibil¬ 


ity of broken U( 1) symmetry forces bj to be expanded 
around its mean-field value = (bj) (which we take to 
be site independent and can always be chosen real) by 
bj = (p + 5bj. If we concentrate on the site at the origin 
6 0 , the Hamiltonian can be rewritten as 


H = H 0 + tf ext + AH, 

H 0 = kn 0 (n 0 - 1) - nn 0 - zJ<j> (b Q + b f 0 ), (4) 

A H = -j£ (Sb\8b 0 + SblSb^ , 

(i,o > 


where i^ ext contains all terms of the Hamiltonian 0 not 
containing the origin “o”. The notation (i, o) means that 
we sum over the nearest neighbors of o, and z is the coor¬ 
dination number. We wish to separate the full partition 
function Z = tr [e _/3// ] as Z = Z ext Z 0 . Here, Z ext is the 
full (and unknown) partition of the system determined 
by the terms in the Hamiltonian not involving the site 
o. It is treated as an irrelevant number in the rest of the 
paper. The partition function Z 0 contains the full local 
Hamiltonian H 0 as well as the correlations introduced by 
“ext” on the origin as follows, 


Z 0 = tr 


e 


p(H 0 + (AH) He J 


(5) 


We approximate the expectation value (A H)n ext by the 
cumulant expansion 


(AH) 


H e: 


-J ('£8*16*0 


(ho) 


H e: 


-\ j2 (E 6 *\ 6 *° E 

< j,o ) 

= 0 - ^6*1 A8* 0 , 


Ha- 


(6) 


where (A H) and (AHAH) are rewritten in terms of 
Nambu operators and A is an unknown 2x2 real-valued 
matrix with two independent components An = A 22 and 
A 12 = A 21 which describes a correction to the common 
mean-field impurity Hamiltonian. The anomalous term 
A 12 , containing processes of the type 5b 2 , is explicitly 
taken to be finite in this notation, since it is known from 
the Bogoliubov theory that deep in the superfluid phase 
it becomes equally important to the normal (diagonal) 
term An, containing the Sb^Sb terms. By ^ we arrive 
at the effective impurity Hamiltonian 


H e = AS* 0 - zJ$ (b a + b\) 

+ kn 0 (n 0 - 1 ) - (m 0 . ( 7 ) 

As can be seen in the Appendix, this effective impurity 
Hamiltonian is equivalent to B-DMFT in the limit 


A(ti - t 2 ) -s- A5(n - t 2 ). (8 ) 

Since A is independent on imaginary time, the Dyson 
equation on the impurity has to be evaluated only for 
a single Matsubara value. The Green’s function which 
mirrors the symmetry relations assumed for A is the one 
evaluated at uj n = 0. The central characteristic of B+U 
theory is that we demand that the condensate and the 
(full) Green’s function of the Bose-Hubbard model eval¬ 
uated at o for zero (Matsubara) frequency coincide with 
the one of system ([ 7 ]) , i.e., 

= (bo) h e 5 ( 9 ) 

g°>°K = o) = - <(*o*S)K = o)) ffB + #t. (io) 
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The paradoxical compatibility of Eqs.(10) and (J8| is spe¬ 
cific for bosonic systems [see also below Eq. (11)]. The 
equations constitute a self-consistency problem, whose 
solution also fixes the factors An and A 12 . This can 
be solved in a unique way if A[G c ,0] is invertible, or, 
technically speaking, if the Luttinger-Ward functional is 
unique. Since the static mean-field (i.e., the decoupling 
approximation with An = A 12 = 0) is always a solution, 
it is easy to convince oneself that multiple (local) minima 
occur (cf. Ref. [26] for a recent discussion). Neverthe¬ 
less, we have been able to determine the physically cor¬ 
rect solution without problem in all parameter regimes 
(see below). In practice, one uses an iteration scheme to 
solve the self-consistency problem. The factors An and 
A 12 follow from the Dyson equation on the impurity site 
at zero Matsubara frequency, 


self-consistency condition ©• <p is thus kept constant 
during the variational calculation of the self-energy. We 
therefore have to minimze 

sn _ sn 5 An sn 5 a 12 

“ <5Au <5X1 + <5Ai2 <5£ ' ^ ' 

We are able to find an analytic expression of since 
the grand potential is defined as U(A, </>) = —lnZ(A,0) 
with Z( A,0) = Tr[e~P HE ^^], giving us 

= <2# 0 - n o ) ff£(A , 0) - |^| 2 , (15) 

—= (2<j>b 0 - b 2 Q ) HE ( a,0) - \4>\ 2 ■ (16) 

After integration, this gives us the relation 


A = SeK = 0) + G C) eK = 0)- 1 -mI, (11) 


H(S, </>) = A( A(E), 0) + B{ A(S), </>) + C. (17) 


with an unknown self-energy He- The connected Green’s 
function on the impurity site G c e is calculated through 
the full Green’s function on the lattice by 

G c ,E(^n = 0) = j^yi J d d kG(u) n = 0, fc), (12) 

and the Dyson equation of the full lattice 

GK = 0 , k)- 1 = GoK = 0 , k)- 1 -^ eK = 0), (13) 

with the bare Green’s function given by Go(cj n , k)~ x = 
[fib — e(fc)] I + icj n cr z , where e(k) is the dispersion relation 
of the lattice and uj n — ^-n are the Matsubara frequen¬ 
cies. 

The approximation (|8| shows that the B+U theory has 
the same functional form as the decoupling approxima¬ 
tion in the non broken phase. In that case only An is 
present but it acts as a shift in chemical potential. For 
the broken phase, An and (i are combined with differ¬ 
ent operators. They control the density of the condensed 
and depleted atoms, whereas A 12 mainly determines the 
anomalous density. According to the Bogoliubov theory 
of the weakly interacting Bose gas, the anomalous prop¬ 
agator is equally important (but opposite in sign) as the 
normal propagator deep in the superfluid phase. In this 
way, the deep superfluid regime is taken care of appropri¬ 
ately in our formalism. The Mott localization is enabled 
by the exact treatment of the density fluctuations on the 
impurity. 


III. VARIATION OF THE SELF-ENERGY 


In order to solve the impurity problem, we consider 
the minimization of its grand potential fi[X,0] with 
respect to its self-energy E, as is also done in self¬ 
energy-functional theory [27] and VC A [16.. The min¬ 
imum with respect to the kinetic condensate term zJp , 
= 0, is already taken care of by the 


with some irrelevant integration constant C and 

^(A(S), 4>) = ^A^An(S), (18) 

0 An 

B{ A(£), 0) = ^Aa 12(S ). (19) 

In order to avoid unphysical results, we have to introduce 
upper bounds on |A^|. From (|6| we see that A cannot 
exceed the kinetic energy of a double hopping process of 
depleted particles, 


^12 


<|A n | <(zJ) 2 (SblSb 0 ). 


( 20 ) 


Furthermore, we require that for all momenta k 
Gn(uj n = 0, &) _1 < —e and det[G(cj n = 0,/c) -1 ] > e 
[where a small e is introduced for stability requirements 
when inverting the 2x2 matrix G(cj n = 0,&) -1 in (13)], 
giving us additional bounds on the self-energy 


Fn > fi zJ + e, (21) 

IS 12 I < V(£n-*J-/i) 2 -e 2 . (22) 


IV. FULL SCHEME AND OBSERVABLES 


By combining Secs. [II] and III we can write down the 
full self-consistent scheme for the B+U theory. Starting 
from an initial guess for <p and A (usually the converged 
mean-field values for A = 0), we calculate a new value 
for p through Eqs. 0 and 0 - Then we search for the 
optimal value of the self-energy by minimizing ( p| ) while 
keeping <p constant. This is done by varying E within 
the bounds (21) and (22) and calculating A(E) by Eqs. 
<0 -(p|), keeping in mind the bound on A ( |2Q| ). Once 
the optimal value E opt is found, the new value for A, 
A(5] op t), is plugged into (m), from which a new value for 
p is calculated. This procedure is repeated until conver¬ 
gence is reached. In the B+U self-consistency all bonds 
adjacent to o are included in He , whereas when comput¬ 
ing the quantities per site, all bonds have to be counted 
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only once. In order to calculate the correct thermody¬ 
namic quantities per site once convergence is reached, 
one therefore has to divide A by 2 , giving us e.g. for the 
density per site (n) = N/V, 


(n) -+ (n o )tf s (A/ 2 , 0 )- (23) 

We can further divide the Hamiltonian into a kinetic [up¬ 
per line in and a potential term [lower line in ([7])], 
giving us expressions for the kinetic and potential energy 
per site 

E kin = -I (A u (5bl5b 0 ) + A 12 (6b 2 0 )) - zJ\cfi \ 2 , (24) 
E po t = ^ (( nl) ~ K)) - (25) 

where the total energy per site is given by E tot = Ekin + 
E pot . It should further be noted that even though we 
do not need to calculate G(r) explicitly in the solver 
and we do not include any retardation in our formalism, 
we can still calculate correlation functions of the kind 
(^(r)B(O)) by 


(A(t)B(0))h e = 7 ?tr 


e -{P-r)H E Ae ~rH E B 


(26) 


or directly in energy space through the eigenvalues of He • 
Since the B+U solver consists of a single impurity, in 
order to compute momentum-dependent quantities, one 
has to resort to the approximate expression 

G(jw„, fc ) -1 = G 0 (*w„, k)~ l - S(*w„), (27) 


with 



FIG. 1. (Color online) (a),(b) Zero-temperature phase dia¬ 
gram of a (a) 3D cubic and (b) 2D square lattice in the vicin¬ 
ity of the (n) ~ 1 Mott lobe calculated with B+U (red dots) 
compared with mean-field (gray line) and QMC [28] (white 
boxes) results. For simplicity, the B-DMFT results are not 
shown here, since they overlap with the QMC data within 
1% M- ( c ) Temperature-dependent phase diagram of a 3D 
cubic lattice with chemical potential ii/U = 0.4 calculated 
with B+U (red dots), compared with mean-field (gray line), 
QMC [28] (white boxes, mostly overlapping with B-DMFT) 
and B-DMFT fTOj (black triangles) results. The B+U results 
for a semi-circular density of states (Bethe lattice, z = 6) are 
shown as a black dashed line. The systematic error bar is 
smaller than the size of the dots. 


X)(ia; n ) — iuj n a z + /il + A — G c e (28) 

By a Fourier transformation this enables us to compute 
such quantities as the momentum-dependent density 

n(k) = -G(k,r = 0 +) - 1 , (29) 


or the critical quasi particle and quasi hole energies at 
zero momentum e p / h which can be evaluated from the 
asymptotic behavior of G(k = 0 ,r) at zero temperature 
through j28f29] 


G(k = 0 , r 



Zp^ T 

Z m e- £ - T 


r -+ +oo, 
r -+ —oo, 


where Z p = Z rn — 1. 


(30) 


V. SIMPLE LIMITS 


theory m , where the self-energy is frequency indepen¬ 
dent as is the case for B+U in our approach. Another 
simple limit of B+U is the Bethe lattice for a semicircular 
density of states given by 

D (e) = |e| < 2ViJ, (31) 

as was also implemented for B-DMFT [10111] , which re¬ 
duces the self-consistency of B+U to one single equation, 

A = -zJ 2 G c ^(uj = 0). (32) 

As can be seen in Sec. ED for the 3D case this leads 
to good agreement with the full self-consistency of a cu¬ 
bic lattice with a much lower numerical cost, since the 
minimization routine described in Sec. m is no longer 
necessary. 


From the relation (20) it is clear that A -+ 0 as J -+ 0. 
Furthermore, also, as U goes to zero, A vanishes, since 
(■ Sbl5b 0 ) -+ 0. Therefore in both cases the mean-field 
limit is recovered. In the case of U <C J the mean-field 
limit is consistent with the weakly interacting Bose gas 


VI. RESULTS 

In Fig. [l] the phase diagram at zero temperature is 
shown for both a three-dimensional (3D) cubic and a two- 
dimensional ( 2 D) square lattice. We compare the results 
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FIG. 2. (Color online) Temperature dependence of (a) the 
superfluid order parameter </> , (b) the kinetic energy F+in, 
and (c) the total energy E t ot for i±/U = 0.4 ((n) ~ 1) and 
U/J — 20 in a 3D cubic lattice calculated with B+U (red 
dots) compared with mean-field (gray line) and QMC [28] 
(black boxes) results. The B+U results for a semicircular 
density of states (Bethe lattice, z = 6) are shown as a black 
dashed line. The systematic error bar is smaller than the size 
of the dots. 



0 Q_._.__ -2 l _,_,_,_,_, 

-ff ^5 5 US' 10 ’ 0 0.2 0.4 0.6 0.8 1 

T/J x/p 


FIG. 3. (Color online) (a) Temperature dependence of the 
local density per site (n) for ii/U = 0.4 and U/J = 20 in a 3D 
cubic lattice calculated with B+U (red dots) compared with 
mean-field (gray line), Bethe lattice (dashed black line), and 
QMC [28] (black boxes) results, (b) Imaginary time depen¬ 
dence of the components of the Green’s function on the impu¬ 
rity G(r) in the superfluid phase for the same parameters and 
T/J m 1 . The normal component Gn(r) = — (6 0 (r)frJ(0)) 
is shown in black dots , while the anomalous component 
Gi 2 (t) = —{b 0 {r)b 0 { 0)} is shown in red squares. 


with mean-field theory, path integral Monte Carlo simu¬ 
lations with worm-type updates (QMC) from Ref. |28j, 
and B-DMFT results from Ref. TO] . The results are 
identical with the B-DMFT results and agree within a 
percent with the QMC data both for the 3D and the 2D 
cases. The results for a Bethe lattice with coordination 
number z = 6 are shown as black dashed lines. As can 
be seen, for the 3D case the simplified self-consistency 
for the Bethe lattice works very well, showing deviations 
only near the tip of the Mott lobe. In Fig. 0c) the 
temperature-dependent phase diagram for jji/U = 0.4 is 



°R X r X M R °0 0.05 0.1 0.15 

(k x ,k y ,k z ) zJ/U 


FIG. 4. (Color online) (a) Momentum dependence of the 
density n(k) at zero temperature for fi/U = 0.5 and different 
values of J/U in the Mott (black, blue) and superfluid phase 
(red, green), (b) Quasi particle and quasi hole energies at 
zero momentum c p /h in the Mott phase at zero temperature 
for B+U (solid lines), QMC 29. (squares), and the analytic 
zeroth-order solution for U —> oo (dashed) for fi/U = 0.5 and 
zero temperature. 


shown and compared to B-DMFT, QMC, and mean-field 
results for a 3D cubic lattice. In this case the lack of 
retardation in the B+U formalism leads to a bigger de¬ 
viation from the B-DMFT results. However, the B+U 
results are still far more precise than the ones obtained 
in static mean field theory. In Fig. [2] the temperature de¬ 
pendence of the superfluid order parameter the kinetic 
energy E^in, and the total energy E t ot is shown and com¬ 
pared to mean-field and QMC for fi/U = 0.4 ((n) ~ 1) 
and U/J = 20. It should be noted that, since the opti¬ 
mization in T, is very sensitive close to the phase tran¬ 
sition, we cannot and wish not to make any statements 
with respect to the order of the phase transition in Fig. [2j 
A local theory such as B+U should be judged for its ac¬ 
curacy on local observables and is by construction unable 
to capture long wavelength physics. Information on criti¬ 
cal phenomena is hence outside its realm of applicability. 
The kinetic energy is very accurate for low temperatures, 
but in the normal phase we find a plateau just as in the 
decoupling approximation. The corresponding local den¬ 
sity (n) per site is shown in Fig. [3] for the same param¬ 
eters, where also the full local Green’s function on the 
impurity in imaginary time G(r) = — (^ > o (r)^>J(0)) j f/ £; 
calculated by (26) is shown for the same parameters and 
temperature T/J = 1. In Fig. pla) we plot the density in 
momentum-space n{k ) = n^(h) + no£/c,o for t/U = 0.5 
and different values of J/U at zero temperature, where 
U'dp(^) are the depleted particles calculated from the con¬ 
nected Green’s function and no = \(j)\ 2 is the condensate 
fraction. In Fig. p^b) the quasi-particle and quasi-hole 
energies in the Mott phase at zero momentum are ex¬ 
trapolated from the imaginary-time dependence of the 
zero-momentum Green’s function through (30) and com¬ 
pared to QMC results from Ref. [29] and the analytic 
zeroth-order solution for U -+ oo. We see both good 
agreement with QMC for finite J/U as with the strongly 
interacting limit for J/U <C 1 
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VII. CONCLUSION 


We have presented the B+U framework for equilib¬ 
rium studies of the Bose-Hubbard model. It captures 
the low-energy physics of a condensed phase as well as 
a Mott localization transition when density fluctuations 
are strongly suppressed on a lattice. The thermodynam¬ 
ics of local quantities of the Bose-Hubbard model can be 
accurately reproduced everywhere in parameter space by 
just three parameters </>, Hu, and £ 12 . By treating the 
self-energy as a variational parameter which minimizes 
the grand potential, B+U can reproduce both the 3D 
and 2D phase transition from the Mott to the superfluid 
phase at zero temperature with an accuracy of about 1% 
near the tip of the lobe and better elsewhere. The B+U 
method can also be applied to finite-temperature systems 


showing a strong improvement on the simple mean-field 
limit but it is less accurate than B-DMFT in locating 
the phase transition line. Just as B-DMFT, being a lo¬ 
cal theory, it can of course never capture hydrodynamics. 
Due to its simplicity and low computational cost, B+U is 
a powerful tool which in the future could be extended to 
clusters in order to study inhomogeneous or topologically 
non trivial bosonic systems, combining exact interactions 
on clusters, embedded self-consistently in a lattice with 
any dispersion. 
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Appendix: Relation between B+U and B-DMFT 


We start from the Hamiltonian |T]), but instead of Z = tr [e @ H ] we treat the partition function in the presence of 
retardation as Z = tr [e~ s ] with the full action S given by 


s = L dT (? b ^ 9 rbi ~ J £ b i ( r ) 6 i( r )j + J o dr Yl 


~^ n i( T )( n i( T ) ~ 1) - M n i( T ) 


(A.l) 


As in Sec. [IT] we expand bj (r) = </> + Sbj(r) around its site- and imaginary-time-independent mean-field value </>, 
giving us 

S = S 0 + S ext + AS, 

So = J dr ^ bl(T)d T b 0 (r ) + ^n 0 (r)(n 0 (r) - 1) - ^n 0 {r) - z J0+(t) + M r ))^ , (A.2) 

J J dr ^2 (sb\(T)Sb 0 (T) + , 


A S = 


and expand the full partition function as Z = Z cxl Z Q by 


Z n = tr 


-s 0 -{AS), 


(A.3) 


As with the Hamiltonian in Sec. [TI]we approximate the expectation value (AS)s ext by the cumulant expansion, 
<AS) Sext » - J^drJ {TS^(t)5^ 0 (t)\ - 1 drdr'J 2 (Y +(r)^ 0 (r) ]T +(t')M>o(t')\ 


(i,o) 


0» 


rP i 

= 0 — / drdr'-5^ 0 {t) A(r — t , ) 5 <£ 0 ( t / ), 

J o ^ 

leading to the effective B-DMFT impurity action PUT] 


(A.4) 


Se = J dr (bl(r)d T b 0 (r) - zJ(f>(b 0 (r) + +r))) ~ j A(r - t')S$ 0 (t') 

+ J dr |yn 0 (r)(n 0 (r) - 1) - Mn 0 (r)^ . (A.5) 
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By the relation 

S[b, 6 f ] = fl dr (6 t (r)a T 6(r) + JJ[6(r), 6 t (r)]), one can see that the action |a. 5| is equivalent to the effective Hamil¬ 
tonian 0 for the ansatz 


A(t - t') -t A<5(r - r') 


In fact, (A.6) reduces the Dyson equation of B-DMFT on the impurity to 

iw n a z + /il + A - T,(iuj n ) = G c _1 (icj r7 


(A.6) 

(A.7) 


If we now take the zero Matsubara frequency of this equation we recover the B+U Dyson equation on the impurity 
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